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Abstract. A Costas array of order n is an arrangement of dots and blanks 
into n rows and n columns, with exactly one dot in each row and each column, 
the arrangement satisfying certain specified conditions. A dot occurring in 
such an array is even/even if it occurs in the i-th row and j-th column, where 
i and j are both even integers, and there are similar definitions of odd/odd, 
even/odd and odd/even dots. Two types of Costas arrays, known as Golomb- 
Costas and Welch-Costas arrays, can be defined using finite fields. When q is 
a power of an odd prime, we enumerate the number of even/even odd/odd, 
even/odd and odd/even dots in a Golomb-Costas array. We show that three of 
these numbers are equal and they differ by ±1 from the fourth. For a Welch- 
Costas array of order p — 1, where p is an odd prime, the four numbers above 
are all equal to (p — l)/4 when p = 1 (mod 4), but when p = 3 (mod 4), we 
show that the four numbers are defined in terms of the class number of the 
imaginary quadratic field Q{\/—p), and thus behave in a much less predictable 
manner. 



1. Introduction 

A Castas array of order n is an n x n array of dots and blanks satisfying the two 
conditions: 

• there are n dots in the array, one in each row and in each column; 

• no two of the (2) line segments joining two dots in the array have the same 
length and slope. 

We shall call an n x n array of dots and blanks satisfying the first condition a 
permutation array, since it corresponds to the well known concept of a permutation 
matrix. A Costas array of order n determines a special type of permutation, tt, say, 
of the numbers { 1, 2, . . . , n }. Specifically, if we have a dot in the (i, j) position of 
the array, we set 

7r(i) = j. 

The second of the Costas array conditions then implies that 

7r(i + fc) — 7r(i) 7^ 7r{j + fc) — 7r(j) 
whenever i j and l<i<i + k<n, l<j<j + k<n. 
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Costas arrays were introduced by the electrical engineer John Costas in the 
1960's. They create ideal waveforms for certain sensor applications, reducing am- 
biguity in interpreting radar and sonar returns. See, for example, [6J, pp. 188-190, 
and [2J. For an arbitrary positive integer n, we do not know any systematic ways 
of constructing Costas arrays of order n, and it is an open question whether Costas 
arrays of order n actually exist for all integers n > 1. Furthermore, in those cases 
where Costas arrays of order n are known to exist, they seem to form a very small 
fraction of the total number of permutation arrays of the same order. Thus, while 
there are clearly 26! permutation arrays of order 26, exhaustive computer searching 
has shown that only 56 of them are Costas arrays. We have a very explicit condition 
for a permutation to define a Costas array, but in practice, it is a condition which is 
difficult to exploit in terms of traditional algebraic properties of permutations, and 
no simple criteria have emerged which enable us to distinguish Costas permutations 
from ordinary permutations. 

In the search for Costas arrays, it is reasonable to ask whether, from a given 
Costas array of order n, we can construct a subarray (or superarray) which has the 
Costas property. For example if the Costas array has a dot at the (1,1) position 
(a so-called corner dot), we can then remove the first row and column to obtain a 
Costas array of order n — 1. This procedure, and similar variants, has frequently 
been employed to obtain Costas arrays of smaller size. Somewhat harder to achieve, 
is to add a corner dot to create a Costas array of order n + 1. An example of this 
procedure is the creation of a Costas array of order 31 from one of order 30. 

The motivation of the present paper, which relates to constructing Costas arrays 
from larger such arrays, is described by the following simple lemma. 

Lemma 1. Let C be a Costas array of order n. Let A be the subset of ordered 
pairs where there is a dot in the (i, j) -position of C . Suppose that there is an 

integer r > 1 such that whenever {i,j) is in A and r divides i, then r also divides 
j . Write n = rm + d, where m and d are integers and < d < r. Then we may 
construct a Costas array of order m by placing a dot in the {u, v) position, where 
V is defined as the unique positive integer such that {ru,rv) G A. 

Proof. We let Vt denote the set of integers lying between 1 and n, and Vl' denote 
the subset of integers lying between 1 and m. Let a be the permutation of 
corresponding to C . We define a permutation a' of 17' by setting 

a\u) = ^ 
r 

We claim that a' has the Costas property. For suppose that 

a' (u + w) — (j' (u) ~ a'{v + w) — cr' [v) 

where \ <u<u + w< m, l<v<v + w<m. It follows then by definition that 

a{ru + rw) — a{ru) — a{rv + rw) — cr(rv) 

But r < ru < ru + rw < rm < n, r < rv < rv + rw < rm < n, and since cr has the 
Costas property, it follows that ru — rv. Hence u = v, which implies that a' also 
has the Costas property. □ 

In the simplest case of Lemma [TJ namely when r = 2, we could (in principle) 
construct two smaller Costas arrays, one by retaining the dots whose coordinates 
are both even, and the other by retaining the dots whose coordinates are both 
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odd. However, such constructions are impossible, as the foUowing theorem indicates 
[HE]: 

Theorem 1. In a Castas array of order n > 2, at least one dot has one even and 
one odd coordinate. 

We will provide examples of larger values of r for which the hypothesis of Lemma 
[T] is fulfilled, although we must admit that the corresponding subarray derived is 
already known. 

Having now discussed our general intentions for investigating Costas arrays, we 
proceed to describe and investigate the field-theoretic constructions of a special 
class of these arrays. 



2. Properties of Costas arrays constructed by field-theoretic 

METHODS 

For certain special values of n, finite fields provide a framework to produce Costas 
arrays of order n, and we intend to show in this paper how the arrays so produced 
relate to the phenomena described in Lemmas 1 and 2. We begin by explaining the 
field-theoretic constructions, as presented in [7]. 

Let p be a prime integer and let q ^ p-f , where / is a positive integer. Let Fg 
denote the finite field of size q and let F* denote the multiplicative group of non- 
zero elements in F^. Since F* is a cyclic group, we may take any two generators 
a and f3 of it (we allow the possibility that a = /3) and construct a corresponding 
permutation array of order q — 2 by putting a dot in the position where 
1 !i J 5: 9 ^ 2, whenever 

a' + 13' = 1. 

It is well known that such an array has the Costas property. We shall call an array 
constructed by this procedure a Golomb- Costas array of order q — 2. 

Suppose that the integer / can be factorized as a product f = de, where d and 
e are positive integers. Write Q — p"^. Then Q'^ and Fg is a subfield of index e 
in Fq . We set 

r = . 

0-1 

Let a be a generator of F* . It is straightforward to prove that a* G Fg if and only 
if r divides i. 

Lemma 2. Let q = Q^, where Q — p'^ is a power of a prime p, and let a and (3 he 
generators of¥*. Let 

r — — . 

0-1 

Suppose that for some integers i and j we have 

a'+(3' = 1. 
Then r divides i if and only if r divides j . 

Proof. This follows from our observation above describing when a' is in ¥q and 
also from the obvious fact that, under our given hypothesis on a and f3, a' S ¥q if 
and only if P' g Fq. □ 
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Lemma [2] implies that ii q — p-^ , where / is an integer greater than 1, a Golomb- 
Costas array of order q—2 always satifies the hypotheses of Lemma [1] The Costas 
array constructed by the lemma is simply a Golomb-Costas array of order p'^ — 2, 
where e is any positive integer divisor of /. Thus the subarray is not new. (We 
note that the case when q = A leads to a vacuous substructure.) 

For each odd prime p, there is also a construction, due to L. R. Welch, of a Costas 
array of order p—1. We define it as a permutation on the set 17 = {l,2,...,p— 1} 
as follows. Let a, c be any generator, any element of F*, respectively, and define 

(J = (Tq .c • ^ — ^ ^ 

by a[i) = j, where j is that element of satisfying 

ca* = j (mod p). 

It is straightforward to see that all these CTq.c are different, and hence there are 

such permutations. Furthermore, each permutation aa^c has the Costas property 
and thus we obtain {p — l)(f>{p — 1) corresponding Costas arrays. We will call these 
Welch- Costas arrays. It is not at all easy to see whether any Welch-Costas array 
has the property described in Lemma [U for a particular value of r. 

3. Parity properties of the dots of a Costas array 

We introduce here a simple concept to generalize the idea implicit in Theorem [1] 
We say that a dot occurring in a permutation array is even/even if it occurs in the 
i-th row and j-th column, where i and j are both even integers. We may likewise 
define odd/odd, even/odd and odd/even dots. In this connection, the following 
notation is convenient to use. 

Definition 1. Let P be a permutation array of order n. We let 

#(e,e), #(o,o), #(e,o), #(o,e) 

denote the number of even/even, odd/odd, even/odd and odd/even dots in P re- 
spectively 

Note that we shall show in Lemma [3] that #(e,o) = #(o, e). In the light of the 
notation just introduced, all that Theorem [1] states is that, for any Costas array of 
order greater than 2, #(o, e) > 0. 

All Costas arrays of order n, where n < 26, are known; there are 143, 635 of 
them. For the convenience of the reader, we tabulate statistical information about 
the location of the dots in Costas arrays in terms of the numbers #(e, e), etc. Table 
[1] displays this information for order n < 15, 15 being the highest order we know 
for which Costas arrays with #(o, e) = 1 exist. We are inclined to speculate that, 
as n oo, #(e, o) ^ oo for Costas arrays of order n. 

Some of the data in the table shows a regularity which can be explained in 
theoretical terms. Thus, for example, rotation of a Costas array of order n clockwise 
through 90° produces another, different Costas array, and a dot at position is 
transformed to a dot at position [j, n+l—i). It follows that if n is even, an even/even 
dot is transformed to an even/odd dot, and an odd/ odd dot is transformed to an 
odd/even dot. This explains a symmetry in the data when n is even. 
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Table 1. Parity populations of dots in Costas arrays of order n < 15 
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Lemma 3. Let P he a permutation array of order n. Then we have 

#(e,o) = #(o,e) 

and 

#(o, o) = #(e, e) + 1, if n is odd; 
ij^{o,o) = #(e,e), if n is even. 

Proof. Let A be the subset of ordered pairs (i, j), where there is a dot in the position 
(i, j) of P. Let £ denote the subset of even integers lying between 1 and n, and O 
the subset of odd integers in the same intervaL There are \£\ elements {i,j) in A 
with i € £ and, by considering whether the second component j is even or odd, we 
see that 

1^1 = #(e,e) + #(e,o). 

Equally, there are \S\ elements (fc, I) in A with / e £ and the same argument shows 
that 

1^1 = #(e,e) + #(o,e). 
It follows that #(e, o) = #(o, e). It is also clear that 

|0| =#(o,o) + #(o,e). 

Since \£\ = |0| — 1 if n is odd, and \£\ — \0\ if n is even, the equations relating 
4l={e,e) and ^{o,o) in the two cases are immediate. □ 

We turn now to the question of enumerating the even/even dots in a Golomb- 
Costas array. 

Lemma 4. Let C be a Golomb-Costas array of order q — 2, where q is a power 
of an odd prime. Then the number of even/even dots in C equals the number of 
non-identity elements z in F* such that z and 1 — z are squares. 

Proof. We may suppose that C is defined using generators a and /3 of F* . Suppose 
that C has a dot at an even/even position (i, j). Then we may write i — 2u, j — 2v 
for unique positive integers u and v, and by definition of C, 

It follows that if we set z = a^", z is a square and 1 — z = (3'^^ is also a square. 
Thus an even/even dot in C determines an element z with the stated property. 

Conversely, suppose that w is a non-identity element in F* such that w and l~w 
are both squares. Then we may write w = a^* and \ ~ w ~ /3^* for unique positive 
integers s and t, with 1 <2s,2t < q — 2. Clearly, 

a + p —1, 

and thus C has an even/even dot in position (2s, 2t). □ 

Corollary 1. Let q be a power of an odd prime. Then all Golomb-Costas arrays 
of order q — 2 have equal numbers of even/even dots, and likewise equal numbers of 
odd/ odd dots, even/ odd dots and odd/even dots. 

We note that there are examples of Costas arrays of order q—2, where q is a power 
of a prime, which do not exhibit the same regularity with respect to even/even dots, 
etc. These arrays, of course, are not constructed by field-theoretic means. 

To enumerate the elements with the property described in Lemma HI let 
denote the subset of F, obtained by deleting and 1. Let S denote the subset of 
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squares in Fj and NS the subset of non-squares. It is well known that lA^^I = 
{q - l)/2 and thus IS*! = {q - 3)/2. 

We now partition into four subsets A, B, C and D defined by 

A = {ze¥l: ze S and 1- ze S} 

B = { z e Fj : z G Af5 and 1 - z e TVS' } 

C = { z e Ft ; z e 5 and 1 - z e A^S* } 

D = { z e Fj : z e iV5 and 1 - z e 5 }. 

Next, we introduce two permutations a and r of Fj defined by 

cr(z) = 1 — z, t(z) = z^^ . 

Provided (? > 5, the two permutations generate a subgroup G of the group of 
permutations of F^, isomorphic to the symmetric group ^3. G permutes the subsets 
A, B, C and D introduced above, and this fact enables us to calculate the size of 
each subset. 

Lemma 5. With the notation introduced above, the following hold. 

(a) If q= 1 (mod 4), then 

<t{A) = t{A) = A, a{C) = L>, t{B) = D. 
Thus \B\ — \C\ — \D\ in this case. 

(b) If q = i (mod 4) , then 

t{A) = C, cr(C) = D, a{B) ^ t{B) ^ B. 
Thus \A\ — \C\ — \D\ in this case. 

Proof. We omit the details of the proof, but note that the key facts are that —1 
is in S if and only if q = 1 (mod 4) and the product of two elements of N S is in 
S. □ 

Since 

\A\ + \C\ = \S\^^^, \B\ + \D\^\NS\^^^, 

Lemma [5] enables us to find the size of each of the four subsets. We summarize the 
details below. 

|A| = (g-5)/4, |B| = |C| = |Z?| = (g-l)/4 

when q = 1 (mod 4) , and 

\A\ = \C\^\D\ = {q~?,)/A, \B\ = {q+\)/A 

when (7 = 3 (mod 4). 

Given Lemma 21 we readily deduce the following information about the number 
of dots with a given parity in a Golomb-Costas array. 

Theorem 2. Let C be a Golomb-Costas array of order q — 2, where q is a power 
of an odd prime. Suppose that the number of even/even, odd/odd, even/odd and 
odd/even dots in C are #(e,e), ^{0,0), #(e, o) and #(0, e), respectively. Then we 
have 

#(e, e) = (9 - 5)/4, #(0, a) = #(e, o) = #(0, e) = [q - l)/4 
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when q=l (mod 4), and 

#(e, e) = #(e, o) = #(o, e) = (q - 3)/4, #(o, o) = (g + l)/4 
when q = 3 (mod 4). 

Theorem [2] shows that, in a Golomb-Costas array, the dots are as uniformly 
distributed into the four positions as it is possible to achieve. Again, Costas arrays 
which are not constructed by such algebraic methods need not show such uniformity. 

It is of interest to enquire whether any such regularity can be detected in Golomb- 
Costas arrays of even order q ~ 2 (in other words, in those defined using a finite 
field of order a power of 2). It is straightforward to show that when g is a power 
of 2, the difference #(e,e) — #(o, e) is an odd integer and examination of all the 
Golomb-Costas arrays of order 2^ — 2 and 2^" — 2, respectively, reveals that this 
integer takes all odd values lying between 35 and —35 in the first case, and all odd 
values lying between 79 and —79 in the second. We have not found any theoretical 
explanation of this numerical data. 

Next, we consider the number of even/even and odd/even dots in a Welch-Costas 
array of order p—1, where p is an odd prime. Here, the results turn out to be much 
more interesting when p = 3 mod 4. Recall that given a generator a and arbitrary 
element c in F*, we have defined a Welch-Costas array of order p — 1 in which an 
even/even dot (2i,2j) corresponds to the equation 

ca^* = 2j mod p. 

This means that if c is a square, the integer 2j is an even quadratic residue mod p, 
whereas if c is a non-square, 2j is an even non-quadratic residue mod p. Similarly, 
an odd/even position (2i + 1, 2j) corresponds to the equation 

ca^^+i = 2j (mod p), 

which implies that 2j is an even non-quadratic residue mod p if c is a square, 
and an even quadratic residue mod p if c is a non-square. It follows that, for a 
Welch-Costas array, we have 

#(e,e)-#(o,e) = ±(i?-F), 

where E is the number of even elements in Q which are quadratic residues mod 
p and F is the number of even elements in fl which are non-quadratic residues 
modulo p. 

Lemma 6. Suppose that p = 1 (mod 4). Then the number of even quadratic 
residues mod p equals the number of even non-quadratic residues modulo p. 

Proof. We use the fact that, under the given hypothesis, —1 is a square mod p. Let 
2j be an even integer in f2 which is a square modulo p. Then —2j is also a square 
mod p and hence p — 2j is an odd quadratic residue in fl. This association shows 
that there are equal numbers of odd and even quadratic residues in fl. Now let a be 
a generator of F* and consider the Welch-Costas array of order p—1 defined by the 
permutation cr of il where a{i) = a'' = j (mod p). We know that #(o, e) equals the 
number of even elements in which are non-squares, and it is straightforward to 
see that #(e, o) equals the number of odd elements in which are squares. Since 
#(o, e) = #(e, o), it follows from our earlier deduction that the number of even 
elements in n which are squares mod p equals the number of even elements in 
which are non-squares mod p, as claimed. □ 
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We now deduce the following analogue of Theorem [2] for Welch-Costas arrays of 
order p — 1 when p = I (mod 4). Note that we have complete uniformity in the 
distribution of dots. 

Theorem 3. Given a Welch-Costas array of order p — 1, where p is a prime with 
p = 1 (mod 4), we have 

#(e, e) = #(o, a) = #(e, o) = #(o, e) = {p - l)/4. 

The main result of this paper is that, for primes p satisfying p = 3 (mod 4), the 
situation is much more complicated than that described in Theorems [2] and [3l We 
first prove a subsidiary lemma. 

Lemma 7. Let p be an odd prime and let V , N be the number of integers in the 
open interval (0,p/2) which are squares, non-squares, respectively, modulo p. Then 
if p = 7 (mod 8), V = E and N — F, where E and F are described before the 
statement of Lemma\d[ If p = 3 (mod 8), V = F and N ^ E. 

Proof. Suppose first that p= 7 (mod 8). Then 2 is a square mod p. It follows that 
if X is an integer in the open interval (0,p/2) which is a square mod p, 2x is an even 
integer in 17 which is a square mod p. Furthermore, any even integer in Q which is 
a square mod p clearly arises by this doubling process. This proves that V = E, 
and the proof that = is identical. On the other hand, if p = 3 (mod 8), 2 is 
a non-square mod p, and then the doubling process maps the integers in the open 
interval (0,p/2) which are squares mod p onto the even integers in fl which are 
non-squares mod p. This proves that V = F in this case, and likewise, N = E. □ 

To state the next result, we need to recall the (complicated) concept of the class 
number of an algebraic number field: 

Definition 2. The class number of an algebraic number field equals the number 
of ideal classes in the underlying ring of algebraic integers. Equivalently, it is the 
order of the ideal class group of the number field. 

For example, the class number is 1 if and only if the ring of integers admits 
unique factorization into prime elements. 

We now apply Dirichlet's class number formula ([1], Theorem 4, p. 346) to deduce 
the following result. 

Theorem 4. Let p be a prime and let h(—p) denote the class number of the imagi- 
nary quadratic field Q{\/~p)- In a Welch-Costas array of order p—1, the following 
hold. 

(a) Ifp = 7 (mod 8), then 

#{e,e)-#{o,e) = ±h{-p). 

(b) Ifp=3 (mod 8), then 

#(e,e)-#(o,e) = ±3/i(-p). 

The class number h{—p) has been much studied and has been the subject of many 
deep results. By a result of Oesterle, h{—p) > 1/55 logp, [9j, p. 361. It is known 
that the largest primes p = 3 (mod 4) for which h{—p) = 1 and 3, respectively, are 
163 and 907, and that there is no prime p = 3 (mod 4) with h{—p) = 2. However, 
the exact value of h{—p) seems to be very unpredictable in general. 
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